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Abstract 

The solution of a radial Schrodinger equation for ^(r) containing a nonlocal potential of the 
form J K{r,r') 'ip{r') dr' is obtained to high accuracy by means of two methods. An application 
to the Perey-Buck nonlocality is presented, without using a local equivalent representation. The 
first method consists in expanding tp in a. set of Chebyshev polynomials, and solving the matrix 
equation for the expansion coefficients numerically. An accuracy of between 1 : 10"^ to 1 : 10"^^ 
is obtained, depending on the number of polynomials employed. The second method consists in 
expanding tp into a set of Ns Sturmian functions of positive energy, supplemented by an iteration 
procedure. For Ns = 15 an accuracy of 1 : 10~^ is obtained without iterations. After one iteration 
the accuracy is increased to 1 : 10^^. Both methods are applicable to a general nonlocality K. The 
spectral method is less complex (requires less computing time) than the Sturmian method, but the 
latter can be very useful for certain applications. 
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I. INTRODUCTION 



The solution of the Schrodinger equation in the presence of a nonlocal potential has 



been the subject of many investigations since 1934 In the optical model one source 
of nonlocality occurs in order to describe the knock-on scattering (a manifestation of the 
identity of nucleons and of the Pauli exclusion principle). Another nonlocality occurs in order 
to describe the dynamic polarization of the target or projectile during the scattering process 
(The Feshbach potential). An interesting study of the relation between both nonlocalities 
and the microscopic structure of target nuclei has recently been presented |2|. The knock-on 
process leads to a semi-separable rank one nonlocality (the meaning of "rank" is explained in 
d Appendix A), while the dynamic polarization nonlocality leads to a general kernel K{r^ r ') 
that acts on the wave function ip in the form of an integral in the Schrodinger equation 

[- + V{r) - k^] i){r ') = - J ^{r, r ')^{r)d^f '. (1) 

The potential V and the integral over the kernel K are in units of inverse length squared, and 
are obtained by transforming them from their energy units into /m~^ units by multiplication 
by the well known factor 2m/^^. Here m is the reduced mass of the incident particle, h is 
Plank's constant divided by 2tt, V is the local part of the potential including the spin orbit 
interaction, and k is the wave number of the incident projectile, related to the center of 
mass energy E by {2m/fi^)E = k"^. For application to the case of nucleon scattering from a 



nucleus, {f? /2m) can be set approximately ec 
of nonlocalities have also been introduced 



described in Ref. {4]. 

For the case that is a semi-separable rank one operator, the solution of Eq. ([T]) was 
initially obtained by various laborious combinations of solutions of a local equation j^, 
subsequently a perturbative method using the Singular Value Decomposition (SVD) [gI was 
developed [3], and later a solution using a spectral ^] expansion into Chebyshev polynomials 
9] was presented [l^. The latter did not make use of perturbation theory but applied only 
to this particular semi-separable rank one nonlocality. Other methods of solution for the 



ual to 20.4 MeV fm ^. Other analytical forms 
3|, and a particular velocity dependent form is 



nonlocality of the type of Eq. ([T]) have also been developed. In 1990 Kim and Udagawa 11 | 



presented an efficient solution for a general nonlocal potential K using a Lanczos iterative 
method. Other iterative methods were also developed jl2|, |l^. Expansions into powers 



of the momentum operator have been presented and local equivalent potentials 



have been developed 



16| and used widely [17|. A different method of solution, employing 



expansions into Chebyshev polynomials, has also been described previously 18| for nonlocal 
potentials. 

It is the purpose of the present article to present an alternate method, different from 
either method described above, by using a combination of the Spectral method and the 
SVD method. The method is nonperturbative, does not depend on a choice of an auxiliary 
local potential Uq, can reach an accuracy better than 1 : 10^^^, is applicable to a general 
kernel K, and its complexity is less than that of Ref. [l8|. A numerical example is given for 
the case of the Perey-Buck nonlocality 16|, and a comparison with an iterative expansion 
into Sturmian functions [l^ is also be presented. That comparison not only serves as a 
check on the Spectral method presented here, but also serves to underline the usefulness 

rn 

of Sturmian expansions. As an example, Sturmian functions were found [20] to be very 
instrumental in describing resonances in low energy nucleon-nucleus scattering, as well as in 
implementing the Pauli exclusion principle, both of which are not generally taken account 
of in the conventional numerical solution of coupled channel equations. 

In section II the version of the Perey-Buck nonlocality kernel used for the numerical 
application will be defined, in Section III the spectral method combined with the SVD 
decomposition will be described, in Section IV a sturmian expansion method for the wave 
function solution, together with an iterative correction method, is described, and Section V 
contains the summary and conclusions. 



II. THE PEREY-BUCK NONLOCALITY 



This nonlocality was first introduced by Frahn and Lemmer 21|, and developed further 
by Percy and Buck [16]. The kernel is of the form 

K{f,r-) = Ui^\f+r'\) -±^^e- ^'^'^''^^ (2) 

where /3 is the nonlocality parameter. For the present application a simpler form will be 
adopted, 

K{f,7^) = VpB{r) h{f,7^) (3) 
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FIG. 1: The Perey-Buck nonlocality function /i/,(ri, r2).for ri = 1.0 fm, for various values of the 
partial wave angular momentum L. 




r, (fm) 



FIG. 2: Same as Fig. [T]for ri = 7/m. 



with 



h{f, r' 



(4) 



7r3/2 /33 

The kernel K given by Eq. (|3]) is convenient for the purpose of numerical calculation but 
is no longer symmetric. For the partial wave decomposition, for each value of the angular 
momentum number L the angular momentum projection /^^(r, r') of the operator h is given in 
the Appendix of Ref. 16| in terms of Spherical Bessel functions and of gaussian exponentials 
in r and r'. The behavior of hL{r,r') as a function of r' (indicated as is illustrated in 
Figs. [T]and|2]for two values of r (indicated as ri), respectively, and using for /3 the standard 
value 0.84 fm. 
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FIG. 3: Potentials Vpb and Vws as a function of r. Both potentials are of the Woods Saxon type, 
given by Eq. ([5]). Vpb is the potential that is multiplied into the nonlocal term hii^r, r') according 
to Eq. ([3]) with parameters given in Eq. ([6]). The one denoted as Vws has parameters given by 
(Vo, R, a) = (— 5/m~^, 9/m, 0.5/m), and is used to define the Sturmian functions described further 
below. 



One sees from these figures that the angular momentum dependence of hilr, r') is reduced 
with increasing values of r, and that hi peaks in the vicinity of r. The decrease of K^lr, r') = 
Vpsir) hiir, r') with distance r or r' is assured by the decrease with r of the factor VpB{f)- 
The latter is taken to be of the Woods-Saxon form 

VpB{r) = + exp[(r - i?)/a]), (5) 

using for the parameters the values 

^ = -3.36 R = 3.5 fm, a = 0.6 fm (6) 

This potential is illustrated in Fig. El where it is also compared with a potential Vws which 
will be used to generate a set of Sturmian basis functions, to be used in a subsequent section. 



III. THE SPECTRAL METHOD 



fl. It 



consists 



A version of the spectral method employed here was developed recently 
in dividing the radial interval into partitions, and obtaining two independent solutions of 
the Schrodinger Eq. ([1]) in each partition. These solutions are obtained by transforming 
Eq. ([1]) into an equivalent Lippmann-Schwinger integral equation (L-S), and solving the 



latter by expanding the solution into Chebyshev functions, mapped to the interval [—1, +1]. 
The corresponding discretized matrices are not sparse, but are of small dimension, and the 
two independent functions are very precise (accuracy of 1 : 10^^^). The solution ip in each 
partition is obtained by a linear combination of the two independent functions, and the 
matrix required to obtain the expansion coefficients has a dimension twice as large as the 
number of partitions, but it is sparse. Details are given in Ref . [9j , and a pedagogical version 



is found in Ref. 22|]. For the present application the division of the radial interval into 
partitions is not made, because the effect of the nonlocal potential would extend into more 
than one partition, making the programming more cumbersome. 

In the presence of a nonlocal potential K, the (L-S) equation for the partial wave function 
ip and angular momentum number L = takes the form 

POO 

^(r) = F(r) + / Goir, r') [V{r')5{r' - r") + K{r', r")]iP{r") dr"dr', (7) 
Jo 

where Qo{r,r') is the Green's function given by 

go{r,r') = -^F{r^)xH{ry), (8) 

where (r, r/) = (r<,r>) if r < r' and (r, r/) = (r>,r<) if r > r', where 

F{r) = sin(A;r); H{r) = cos{kr) + ism[kr), (9) 

and where k is the wave number. 

The nonlocal part K of the potential is expanded by means of the Singular Value 
Decomposition (SVD) as follows. First a set of Chebyshev support points ^j, i = 1,2,..N 
are defined in each partition and the corresponding discretized matrix ICij = K{^i,^j) 
is obtained, where N is the number of Chebyshev polynomials Tq, Ti, ...T/v-i to be used 
in the expansion. These support points are the zeros of the first Chebyshev polynomial 
T/v(a;) not used in the expansion, mapped into the radial partition interval, as described 
previously Q. However, in the present apphcafou only o,re large ^artitiorr .s used. 
The singular value decomposition of a x matrix /C is given by [6], /C =USV'^, or 
/Cjj = ^^,s<^s('^^i,s)^ hj = l,2,..,N. This result can be expressed in terms of the 
column vectors Ug and Vg of the N x N unitary matrices U and V, respectively 

N 

/C,, =J2'^s{^)crs[vsm\ (10) 

s=l 



and S is a diagonal matrix containing the singular values as, with s = 1, 2, .., N. The as are 
positive numbers, ordered in descending values, and the symbol f signifies transposition and 
complex conjugation. The matrices U and V"^ are unitary, but are not orthogonal to each 
other. 

According to the Eq. fllOp the matrix JC can be expressed in terms of a sum of products 
(u) X {vY, each of which represents a rank one N x N matrix. Thus, if the values of the wave 
function ip at the support points were expressed as a column vector i/j, with ipj = ip{C,j), 
then the integral in over dr" of the dot-product [vs]'' ■ il^ , as required in Eq. ([7]), would be 
given by 

N 

{Vs\ -0) = (11) 

i=i 

where wj are integration weights. Hence the action oi K on ip can be broken into simpler 
expressions of the type 

N 

(/CV), = 5^ Us{^)as{vs\1^^). (12) 

s=l 

In the spectral method, the integral (Vslil^) is done by mapping ?/'(r) — )■ iplx) onto the 
interval [— 1 < a; < +1], by expanding iplx) into Chebyshev polynomials, 

AT 

^(^) = 2Z (^jTj-ii^)- (13) 
j=i 

and by using the integral properties of the Chebyshev polynomials. The various overlap inte- 
grals involving Chebyshev polynomials are obtained by a method denoted as the Chebyshev- 
Gauss integration procedure, and is given by Eq. (50) of Ref. 22]. By this means Eq. ([7]) 
can be transformed into an equation containing matrices acting on the expansion coefficients 
ttj, and its algebraic solution permits one to obtain the coefficients aj, and hence ?/'(r) via 
Eq. ( JT3l) . The method utilizing the SVD expansion, described above, can be utilized for any 
non-singular matrix, and hence the method applies to any non-singular kernel K. 

The solution ip of Eq. ([7]) for k = 0.5 fm^^ is displayed in Fig. HJ The effect of the 
nonlocality is to "push out" the wave function to larger distances. This can be seen from 
Figs. Oandini which compare the real and imaginary parts of the solution ip for k = 0.5 fm~^ 
in the presence of the nonlocality with the solution for which the non-locality, given by Eq. 
(jl]), is replaced by a delta function, h^^r, rf) 5{r — r'). 
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FIG. 4: The wave function ip{r), solution of the Schroedinger eq. ([7]) with the nonlocal potential 
kernel K given by Eq.([3]). The solution is obtained by the spectral Chebyshev expansion method 
for k = 
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FIG. 5: Comparison of the real parts of the local and nonlocal wave functions ^/J. The nonlocal one 
assumes the presence of the Perey-Buck nonlocality, the local one replaces the nonlocality hL{r,r') 
by 6{r — r'). Here /iL is defined in Eq. ([4]) with L = and k = 0.5 fm^^ 



A. Accuracy analysis 

The error of the spectral method for the nonlocal potential K is obtained by comparing 
solutions for three different values N = 51, 71, and 301 of the number of Chebyshev expansion 
functions used in the radial interval [0,20 fm] with the solution obtained with = 501. 
The results are shown in Fig. [71 and listed in Table HI which also displays the respective 
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FIG. 6: Same as Fig. [5] for the imaginary parts of the wave functions. 
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FIG. 7: Estimate of the error of the spectral method, by comparison of the wave function ip for 
various values N with each other. Here is the number of Chebyshev polynomials used in the 
solution of Eq. ([7j) in the radial interval [0, 20 fm]. Plotted are the absolute values of the difference 
of tp between the results for N = 51, = 71 and N = 301, respectively, with the result for N = 
501. 



computing times. The calculations are done in MATLAB, on a desktop using an Intel TM2 
Quad, with a CPU Q 9950, a frequency of 2.83 GHz, and a RAM of 8 GB. Had the radial 
interval been subdivided into partitions, a further decrease of computing time would have 
been achieved. Such partition division is required if the wave function extends out to large 
distances, as was the case, for example, in the calculation of the Helium di-atom bound state 



231] in the presence of a local He-He potential. 
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N 


accuracy 


time{s) 


51 


2 X 10-^ 


0.4 


71 


10-7 


0.7 


301 


4 X IQ-ii 


13 



TABLE I: Accuracy and computing time of the spectal method. N is the number of Chebyshev 
support points in the radial interval [0,20 fm] 
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FIG. 8: The singular values ag, s = 1,2, ..30, for the Singular Value decomposition of the nonlo- 
cality kernel Ki{r,r') = VpB(r) hL{r,r'), for L = 0. The function hi is defined in Eq. (A. 2) in 
Ref. [16] and is illustrated in Figs. [T]and[21 Vpb is defined by Eqs. and ([U]), and is illustrated 
in Fig. [3j The number of Chebyshev support points, which is equal to the number of Chebyshev 
polynomials used, is = 302. 



An illustration of the singular values cr^ for the case of the Perey-Buck nonlocal kernel, 
using only one partition from r = to Tmax = 15, with = 301, is shown in Fig. [HI One 
sees that for s > 30, the corresponding values of as are less than 10~^, and the expansion 
flTU]) could have been truncated at s = 30 if an accuracy of 1 : 10"^ had been sufficient. The 
Chebyshev expansion coefficients aj, j = 1,2, ..A^ of the wave function ip, Eq.l TT^ mapped 
into the interval [— 1,+1] are illustrated in Fig. [9l The magnitude of these coefficients 
decreases exponentially with the index j, reaching the value ~ 10~^^ for j = 100 and 
beyond. Had the expansion been truncated at j = 30, the error of the wave function would 
have been 1 : 10~^. 

In conclusion, the accuracy of the solution of Eq. ([7]) with the spectral expansion method 
for a nonlocality described by Eq. ([3]), and using A^ = 51 Chebyshev polynomials is 1 : 10~^, 
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FIG. 9: The Chebyshev expansion coefficients aj of the wave function ip that satisfies the nonlocal 
Eq. ([7]), where K{r,r') is the Perey-Buck Kernel ([3]). The radial interval is [0, 15/m], the wave 
number k = 0.5 fm~^, and the total number of Chebyshev polynomials is 302. 

while using = 301 polynomials, the accuracy is better than 1 : 10^^°. 

IV. THE STURMIAN EXPANSION METHOD 

The method consists in expanding the unknown solution ip of Eq. ([7j) into a basis set of 
"global" functions, formulate an equation for the expansion coefficients a^, and if N is 
not large enough, iteratively correct for the truncation error. The expansion into Sturmian 
functions was described in Ref . [l^ , and only a few basic equations relevant for the present 
case will be repeated here, while Appendix B contains further details. References to many 
other applications of Sturmian functions to the solution of physics problems can also be 
found in Ref. Q. 

For the present application the operator O in the general one- dimensional integral equa- 
tion 



JO 

where it is assumed that the local part V of the potential in Eq. ([7]) has been set to zero. 
The function F is the driving term and Qq is the Green's function, described in Eqs. ([9]) 




(14) 



to be solved for ip is, in view of Eq. ([7]), given by 




(15) 



11 



and (fTTj) . respectively, both assumed to be known. The shorthand form of Eq. (fT^ is 

ip = F + Oi). (16) 



The iterative solution of Eq. f lHj) is achieved by approximating the operator (9 by a separable 
representation of rank A^, defining the remainder A^^"* as 

/\n = 0- On, (17) 

and iterating on the remainder. 

The approximate discretization of the kernel O into a representation of rank is ac- 
complished by using a set of auxiliary positive or negative energy sturmian functions $<i(r), 



s = 1, 2, .., iV and is of the form 



N 



0N{ry)=J20'^s)^^^{'^sV. (18) 

Here the symbol ) denotes that the quantity to the left of it is evaluated at position r, and ( 
denotes that the quantity to the right of it is evaluated at r'. The bra-ket {^gV^s) denotes 
the integration {^gV^g') = ^s{T")V{T")^s'{T")dr where ($s is not the complex conjugate 
of $s, and V{r) is the local potential used in the definition of the Sturmians. 

The sturmian functions $s are eigenfunctions of the integral kernel go{r,r')V{r') 

POD 

Vs^sir)= goir,r')Vir')^sir')dr', s = 1,2,3,.... (19) 

^0 

with rjs the eigenvalue, and V{r') the sturmian potential. The differential Schrodinger equa- 
tion corresponding to Eq. ([T^ is 

(dydr^ + E) =A,F$„ (20) 

with As = l/f]s- The Sturmians for positive energies are not square integrable, but they are 
orthogonal to each other with the weight factor V (that is assumed to decrease sufficiently 
fast with r). The normalization of the Sturmians adopted for most of the present discussion 
is 

{<l>sV<l>s')=VsSs,s'- (21) 

Because of the completeness of the sturmian functions, one has the identity 

oo _ 

6{r -r') = J2 '^^(^) (22) 



s=l 



12 



which shows that the larger the number of terms in the expansion ( !T8|) . the better is the 
approximation of On to O, provided that O is compact. 

The first step in the solution of Eq. (ITBl) is to obtain the solution of the approximate 
equation 

/oo 
ON{r,r")F{r") dr". (23) 
n n 

The solution can be obtained algebraically |20[ by making the ansatz [I9j 



N 



s=l 

and the coefficients Cg , 3 = 1,2, ..N , are obtained from the solution of the matrix equation 

^ 1 



where 



Eqs. ( I25ll and ( 1261) are obtained by inserting (1241) into ( l23l) . making use of Eq. (iTSjl . 
multiplying the resulting equation by [l/($sV ('^'s('^) ^^(?") on both sides, integrating 
over r, and making use of the normalization fl2T]) . The matrix element involves 
a triple integral, while the result {(^sK^^s') in Eq. (126!) requires only a double integral in 
view of Eq. flTI?]) . This simplification is one of the advantages of using Sturmian functions 
in the expansion (ITSl) of (9. The integrals in the equations above are carried out using a 



Gauss-Chebyshev procedure |22|] that has high accuracy, and does not require the use of the 
SVD decomposition. 



Had the expansion Yli^=i Cs\Qo^^s) in Eq. (12^ been replaced by the expansion 



X^fLi ds\^s), then the relation between the two sets of coefficients would have been 

N 

ds = Y,Ms,s'Cs' s = l,2,..N. (27) 

s'=l 

Since N can be larger than A^, the expansion in Eq. is preferable because (9 $s)r niay 
lie outside of the space spanned by the functions $s)r, and hence is more general. However, 
for the present application N = N. 

Results for the function J^, given by Eqs. (p^ to (126|) . were compared with the wave 
function obtained by the Spectral Chebyshev method described in section HI. The difference 
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FIG. 10: The error of the sturmian expansion of the solution of Eq. ()23p . for different numbers 
N of Sturmian basis functions. The auxiUary Sturmian potential V, is of the Woods Saxon form, 
given by Eq. ([5]) with parameters (Vo,i?, a) = 5/m~^, 11/m, 0.5/m). The radial interval is < 
r < 20 fm, and the wave number is k = 0.5fm~^ 




Ns 


Np 


T 




10 


301 


2 X 10-^ 


7 X 10-^ 


15 


301 


9 X 10"^ 


2 X 10"*^ 


20 


453 


3 X 10"^ 


3 X 10"'^ 



TABLE II: Accuracy of the Sturmian expansion. Ns is the number of Sturmians, Np is the number 
of support points in the radial interval [0,15 fm] 

between the two results was taken as a measure of the accuracy of the Sturmian expansion, 
that in turn depends on the range R of the Sturmian auxiliary potential V, and the number 
of sturmian functions used. Here R is defined in Eq. (|5]). The dependence of the accuracy 
of J-" on the number A^ of sturmians for a fixed value of i? = 11 fm is illustrated in Fig. [TOl 
It shows that with A^ = 20 an accuracy of 1 : 10~^ is achieved. The accuracy for different 
values of R and A^ is is summarized in Table HIl and is illustrated in Figs. [12] and [131 The 
computational time, carried out in MATLAB, is listed in Table IIIII This result does not 
include the time required to calculate the Sturmian functions. 
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Ns 


Np 


T{s) 


V^iter.(s) 


10 


301 


11 


0.04 


15 


301 


11 


0.04 


20 


453 


26 


0.11 



TABLE III: Computing time of the Sturmian expansion. Ns is the number of Sturmians, Np is 
the number of support points in the radial interval [0,15 fm] 

A. Iterative corrections to T. 

The iterative correction to the Sturmian expansion of J-", Eq. (!24|) is carried out by 
defining the remainder according to Eq. f lT7|) and iterating on the remainder. If the 
norm of A^v is less than unity, the iterations should converge |l9j. Since the numerical 
complexity of performing iterations is less than the complexity of solving a linear equation 
with a matrix of large dimension, this method can be computationally advantageous, as 
is shown in the 4*'^column of Table llllt and furthermore the exact eigenfunctions of the 
operator O need not be known. 

The iterative corrections to ip proceed according to 



where the Xn+i are related to Xn through the iterative equation [19 1 

Xn+l = ONXn+l + ^NXn, n = 0,l,2, .., 

with xo = J^- If the expansion of Xn+i is given by 

N 



(28) 



(29) 



(30) 



s=l 



then the coefficients ds^"'~^^^ obey the algebraic equation 



N 



J2 iSs,s' - Ms,s') dj 



(n+1) 



s'=l 



{^sV AnXu), 



(31) 



where is given by Eq. 0261] . However, before inserting the numerical value of Xn{f) 



into the right hand side of Eq. fl3T]) . Eq. (!30l) is replaced by ci"^ l^o^'^'s), with 



Es'=iiM'')s,s'dj!''\ according to Eq. ([2 
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FIG. 11: Iterative corrections to the wave function ■0, using 10 Sturmians, which are described in 
Fig. [TOl The result labeled F illustrates the value of J-" as described in Fig. [TOl The other results 
are obtained after 1, 2 or 3 iterations. The maxima of the three curves occur at 3 x 10^'^, 10~^, 
and 4 x 10^^, respectively. 



The results for the number of Sturmians = 10 is illustrated in Fig. [TT]for R = 11 fm. 
It is seen that as the number of iterations increases, the result for ip converges but not, 
within the accuracy of the calculation, to the value obtained from the spectral method. 
This is in contrast with the case of a local potential. There the converged result of the 
iterative Sturmian method was found to be in good agreement with the spectral result [3] 
of 1 : 10~®, as compared with the less accurate result for the present nonlocal case. It is 
suspected that the reason for this difference is related to a difference in the treatment of the 
long-range part of the potential, as is discussed further in Appendix B. On the other hand 
the error of (compared to ip) is much smaller in the present nonlocal case than it is for 



the local case 



191]. That difference is very likely due to the presence of a repulsive core in 



the scattering potential for the local case, but is absent in the nonlocal case. 

The effect of the range R of V and the number of Sturmians is displayed in Figs, 
and [131 These figures show that for a given A^, the accuracy of J-" is better for a smaller 
value of R, but the iterations give a more accurate value for ip for a larger value of R. Table 
mil shows that the iteration time is considerably less than the time to compute J-", and in 
order to obtain the same accuracy, the latter is considerably longer than the computing time 
for the Spectral Chebyshev expansion. 
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FIG. 12: Maximum absolute error of tp, as a function of the number of Sturmian basis func- 
tions. The Sturmian potential V is defined in Eq. ([5]) with the parameters (Vo,i?, a) given by 
(—5 fm'^'^, 9 /m, 0.5 fm). The line labed "-F" represents the result for T, Eq. (j24p . while the 
line labeld "iter 1" is obtained by correcting by one iteration, as described in the text. 
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FIG. 13: Same as for Fig. [12] for a larger range R = 11 fm of the sturmian potential. 



V. SUMMARY AND CONCLUSIONS 



The solution of a radial Schrodinger equation for ip{r) containing a general nonlocal 
potential of the form J K[r,r') ip{r') dr' is obtained by means of a spectral expansion 
into Chebyshev polynomials {q], combined with a Singular Value Decomposition. For a 
semi-separable kernel K of rank one, that occurs in exchange scattering due to the Pauli 
exclusion principle, a spectral expansion method has been previously devised [lO|. But for 
a general K the present method is quite different, and has not been presented before. The 
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results for a numerical example of K given by a Perey-Buck ansatz [l^ are calculated and 
their accuracy, as well as the required computer time, is investigated as a function of the 
number of Chebyshev polynomials employed (see Table [T]). For a wave number k = 0.5 fm~^ 
an accuracy of between 1 : 10^^ to 1 : 10^^° is obtained as the number of polynomials in the 
whole radial interval is increased from 50 to 300. A second method is presented, based on 
an expansion into Sturmian functions of positive energy, supplemented by an iteration 
procedure For = 15 an accuracy of 1 : 10~^ is obtained without iterations, and 

after one iteration the accuracy is increased to 1 : 10~^, as shown in Fig. [131 The iterations 
converge quickly, but not to exactly to the value given by the spectral method. It is suspected 
that the reason is due to the range of the Sturmian functions employed not being as large 
as the range of the nonlocal potential K, as is further discussed in Appendix B. 

The method, being applicable to a general nonlocality K, opens the way to formulate 
optical potentials that incorporate the physical effects that are the source of the nonlocality. 
This study is particularly relevant for Astrophysics, where the colliding nuclei are generally 
unstable, and hence the resulting nonlocalities are more pronounced. 

Acknowledgment: The author is grateful to Professor M. Jaghoub for rekindling his in- 
terest into the subject of nonlocalities. 

Appendix A: The rank of a matrix 

Various definitions of the rank of a matrix can be found in Chapter 9, p. 432 of Ref. 



2J]. According to one of the definitions, the rank of a matrix A is the order of the largest 
square sub-matrix of A whose determinant in not equal to zero. A practical definition is 
obtained via the SVD decomposition of the matrix A according to which the rank of A is 
the number of non-zero singular values of A. For example: given two column vectors u and 
V, both having A^ elements Ui and Vi, i = 1,2, ..N, then v'^ -u = X^ili '^i'^i is a number, while 
M f ^ is a A^ X A^ matrix W. Here the superscript T means transposition. The matrix W has 
only one nonvanishing singular value, hence its rank is 1. The 4x4 matrix A given on p. 



433 of Ref. 2J] has 3 non zero singular values, hence its rank is 3, in agreement with Ref. 



24|. 



In the case of exchange scattering, the non locality K is given by K{r, r') = f{r)g{r') 
for r/ < r, and g{r)f{r') for r' > r. Once the radial distance r is discretized into A^ values 
Vi, the function of r- is interpreted as a line vector, while the function of rj is taken as a 
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FIG. 14: The Sturmians eigenvalues of Eq. (j20p obtained with a Woods-Saxon potential Vws 
described in the caption of Fig. \T2\ with a wave number k = 0.5 fm~^ 



column vector, and the kernel K is the separable combination of two rank 1 matrices, each 
of dimension N x N, and hence it is called semi-separable of rank one. 

Appendix B: Sturmian functions 

Some properties of Sturmian functions [19] will be recalled here. According to Eq. ( !T9|) 
Sturmian functions are defined in the radial interval [0,oo]. However, the numerical eval- 
uation of (fT9l) has to be carried out only in the interval < r < -Rmax, beyond which the 
Sturmian potential is negligible. An example of the spectrum of the A eigenvalues for the 
Woods-Saxon potential Vws illustrated in Fig. [3|, for a wave number k = 0.5 fm~^, is shown 
in Fig. [TH The imaginary parts are slightly negative, while the real parts are monotonically 
positive. This is in contrast with the case of a sturmian potential that has a repulsive core, 
for which some of the eigenvalues have a negative real part, and a correspondingly positive 



imaginary part 191]. 



The radial r— dependence of some of the Sturmian functions is illustrated in Figs. [15] 
and [m Their normalizations differ from Eq. (I2T]) so that all acquire the same Hankel 
function asymptotic behavior, since they obey Eq. f|T9|) . The Sturmian potential V has a 
range R = 9 fm, and the figures show that for r > 8 fm, these Sturmian functions are not 
strongly linearly independent from each other, and hence cease to form a practical complete 
expansion set. This is the reason why the Sturmian expansion for a Woods-Saxon potential 
with a range R = 9 fm does not give a result as precise as the expansions for R = 11 fm 
and even that range may not suffice to obtain an accuracy better than 1 : lO'', displayed in 
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FIG. 15: The real parts of Sturmian eigenfuncions of Eq. (|20p obtained with a Woods-Saxon 
potential Vws described in the caption to Fig. [12l whose eigenvalues are illustrated in Fig. [HI 
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FIG. 16: Same as for Fig. [T5l for the imaginary parts of the Sturmian functions. 



Fig. [131 This conclusion is further corroborated by examining the magnitude of the kernel 
K in the region where the Sturmian functions loose most of their independence, which is of 
the order of 10~^ /m~^. 
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